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WORKSHOP PROGRAM -- Multigrid and Molecular Dynamics 
October 10-12, 1995 


Dept. of Applied Mathematics & Computer Science 
The Weizmann Institute of Science, Ziskind Building, 
Lecture Room 1 (basement level) 


Note: The program may change even during the course itself according to 
the interest of participants. 


Tuesday, October 10th 


8:30 - 9:15 Registration, organization, coffee 
9:00 - 9:15 Introductory remarks 
9:15 - 10:30 Introduction: tasks and types of multiscale 
computation (Brandt) 
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11:00 - 12:15 Multiscale decomposition of forces and fast 
force summation (Brandt) 
12:15 - 13:30 Lunch 
13:30 - 15:00 Elements of linear multigrid for PDEs (Stueben) 
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15:30 - 17:00 Poster session, including: 

-Gomathi Ramachandran: Buckling transitions in 
superhelical DNA: dependence on the elastic 
constants and DNA size 

-Bimalendu Mishra: Error analysis in Langevin dynamics 
simulations 

-Margaret Mandziuk: Resonance in the dynamics of 
chemical systems 

-Eric Barth: LIN and a family of related methods for 

molecular dynamics simulation 

~Hongmei Jian: Computer simulation of DNA using a 
discrete worm-like chain model 

-Gerd Winter: BEMOLPA project 

:30 Mini-vans depart for dinner at Pearl of the Sea 
restaurant, in Rishon L’Tzion 


Wednesday, October ilth 


8:45 - 9:00 Coffee 
9:00 - 10:30 Multigrid Monte-Carlo and stochastic coarsening 
(Brandt; Mack?) 
10:30 - 11:00 Coffee 
11:00 - 12:15 Global and discrete-state optimization: multiscale 
annealing (Brandt) 
12:15 - 14:00 Lunch 
14:00 - 15:00 Full multigrid and nonlinear multigrid (Stueben) 
15:00 ~ 15:30 Eric Barth and Margaret Mandziuk: Molecular dynamics 
Simulation at large timesteps 
15:30 - 16:00 Coffee 
16:00 - 16:30 Bob Skeel: Efficient use of fast electrostatics in 
molecular dynamics 
16:30 - 17:00 Chris Lambert: Efficient dense Hessian computation in 
molecular minimization 
17:00 - 17:30 Oren Becker: Reduced variable molecular dynamics 
17:30 - 18:00 Dexuan Xie: A remark on algebraic multigrid analysis 
19:00 Dinner in town at one’s choice 


Thursday, October 12th 


8:45 - 9:00 Coffee 
9:00 - 10:00 Algebraic multigrid (Stueben, Ruge) 
10:00 - 10:30 Coffee 
Multiscale methods in molecular dynamics (Bai, Brandt) 
Discussion on future directions; box lunch 
16:00 Bus departs for dinner in Yaffo, and evening at Tel Aviv 
seaside promenade 
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Introduction: Task and ‘l’ypes 
of Multiscale Computation 


oon particles : nO 10”, Syed 


® Calculating O(n?) forces 
at each time step 


e Extremely short steps ° o—"" 5 


e Wreng attraction basins 


e thermal statistical fluctuations : 
Need represent many probable Uo. 


All addressed by 
multi-scale methods 


multigrid multi— level multi-resolution 
ee Renormalization Group , 
Wavelets, Multipele, Fractals 


Hardware - 
Ever smaller and Faster 


Past: lo times. faster 
each ~6 years 


Future: Pavalle| piconets 


Algorithms 


Required amount of computations 
rise slower with 
increased problem s(ze 


Past: Compara ble contribu bron 


future: Malti scale methods 


Highly compatible 
Exee\ at large problems 


| 
| 

| 
i. 

| 


OO «gt. Ge " of of @& | 
| fot a a Ge GP Be | 
| Ce ae ae | 
| 2 aL ZO ch | 


& 
C 
\s 
( 
\s 
& 


Q, 

2 

‘\ 
. Cc 
Q 


YS 


{> 
| 
= 
—_ 


| 


4 


@ Similarly : For non-uniform atom distribution. 


e Collective force : via the finer-grid forces 


Eneray ( grid) 


Multigrid cycle 


Sweeps on 


scale ga (particles ) 
grid h 
grid 2h 


Very large scale changes in one cycle 


Negligible work On coarse levels 
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FULL MULTIGRID (FMG) ALGORITHM 
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Ne pre 
multigrid eye leV 
} interpolation (order l+p) 
to a new grid. 


interpolation (order m) residual transfer 
of corrections : 
6) dD relaxation sweeps 


enouph sweeps 
O algebraic error Ge) si F 
or direct solver 
< truncation error 


Fully efficient ae minimal work units 


Multig rid Selvers 
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Linear scalar elliptic equations (19%0)* 
Nonlinear 9 Adaptive grids FAS 
General boundaries : Bc ™ 


Discontinuous coefficients vlavu)=f 
Disordered coefficients AMG 

Non scalar PDE systems” elasticity , NS, 
Non- ellipticity convection high Re 
Indefinite : waves <—> rays 
Tnverse problems A=% DA 
Time dependent | space + time 
Dense NN motrices O(n) 

Crucial features invisible to coarse grids 


Gauge freedom | 
Dirac eqs. 


Topolegies 


Near zero modes , 


* Rigo rvs 


Cowpressible Newer - Stohes i ee 


“MAU + QUU,+ GUY - Pray) (Ut) +P, =O 
~Ppbv + guy, + PUY, - +p) (arty )y + Py =O 
-uAE +g & + gu Ey + p(Yx +4) 

pietUy)* ~ a (a+)? -apnlugtyy) =O 


(9 u), + (ev) y ~VASG =O we ES 
P = p (€, 9) ( | Gy) : 
\ 
iy-ellsptic BVP . anes - 
ache a Me ay 
S ia / 
() we ye 


: 4 | . “Sao 

KR VIS Cy CcCase€: a /» “ Pad gt rmax (fal, [4) 

£ = length ak whieh (u,y; €) change . > Euicae 64 
Usv.elly thir, are VISCOUS loyers , 


Tncomressible: Q 0 %\/¢ 
uw} =Q 


o Q Oy 
Ox ‘Oy ) P 
det ie ze Fa) Q Ellvabee | 


Q= Ot 42, +0, . 


h- principal Lust det LL 


Cauchy - ay 2 ie ; 

Riemann & 2t)( ‘ = 0 ras 

20 -A 0 Ox U 

Stokes oO -A- ay if Oo ne 

Ox ey 0) p 

2D 

Incompres SI ble : ; a u —-QA 

Navier - y]}\ v | 

Stokes dy Ry 0 P Qs-SAtUwY 
puvY Oo 9 Qy 

oe 


0 0 ay - R 1 


Compressible Navier- Stokes ku (2m4r)d (u-V) 


(on he viscous scale ) 


Central Cauchy - Riemenn an 
Central (Navier-) Stokes Q" ne 


Po= (Pcr»Péa> Pes) 
P; 


rij= |Pe-P; | 


E°(p) = E Cpyapaoe) = Va (neil) 
E° (p) = a E°(p) 


With external forces: E(p)=E (p)-2 fi-pi 
Find p* such that E(p*) = min E(p) 


Particle - by = particle minimization 


=> Extreme slew down 


Needed: collective mioves 
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(for smooth uh) 


e Minimize E(u), only then displace particles 
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=> Similarly on coarser grids 


=> Elasticity — Plasticity — Strain limits 
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Minimization of E(u) 
Point ye point mintmizatron 


—s slow down ! 
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Ou qrrich problem 
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Fine scale Cearse scale 


Uniformly alliptic eq. pame a: 


\Wave eq. "Rays" eq. 
AU + Kr =O) 

Particles Continuum eqs, 
Lennard ~Jones Elasticity 
ener ou eqs: 


Mn Wat eecl On 


Derivation of macroscopic ‘equations 


(Simpler than senerwmelleacion group ) 


Stochastic Deterministic 


Coe 
Determintstie Stochastic 


AL 


Oet wmizZart IOV min E(u) 
a Fast convergence Neat optimum 


Global optimization | 
Escape false attraction basins 


multi-scale attraction basins 
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min E(e) 


False attraction basins 


Particle by particle minimization 
trapped in local attraction basins. 
~e(u) /T 


Simulated annealing P(u) ee 
TS O 


trapped in large-scale attraction . basins, 


Multi-level annealing 

Collective moves, at ell scales 

conver ges (fast, in probability) . 

o A large-scale move is decided only after 
optimizing around it at all finer scales. 

@ Fast anncaling ot each level. 


2 Recombinations ak each level. 


aH 


Thermal Statiskica| Fluctuations 

4 Need to produce (stochastically, in a 
correct probability distribution ) 
independent configurations 


2. Need to produce very many of thew 
(to avev-age out deviations) 


Multiarid Monte - Carle 

‘a Collective moves, at a\\ scales 

2. Productng many . samples on coarse grids 
(lcktle werk per sample , averaging cub large~ scaly 


fluctuations s loca\ Cluctuakicns ave self ~ averaging 
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Optimal : 
Thermodynamic [ae see 


in Ofe=) computer operations 
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© General non-linear PDE 


steady stote: Elliptic, non- elliptic 

time dependent 

inverse problems... Optimal control 

Waves <—> rays. Many- eigenfunction ; 
topological singularites + disorder (Dirac es.) 


o Integral transform, equations 
® n- body interactions 


@ Mo lecular mechanics - Ground stales 
Equilibrium 
Dynamics 


o Global optimization 
e Linear Programming 
e Monte Carlo 

° Peet aes 

° Multiple - fermion path integrals 

e Derivation of macroscopic equations 
e [mage processing 

e Tomography X rays, MMR, PET,...radar 
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Decomposition of Forces 
and Fast Force Summation 
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Intearal Transforms 
V(x) = Joy U(y) ely 
X= (x, X2,X3) . ‘ie (41,42, Ys) 


Asymptotically smooth kernels 


G(x,y) a e* *d Fourter 
| od Laplace 
ek /o* Gauss 
—ke 
[x-4| ‘ fog [xy Potential 


Oscillatory kernels 
E(x,4) otk Ix-yl 


Integro -diffevential equations 
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G (x,4) = Csi 0 ee G (x,y) 


directly coarser grid 
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Car: (s) (r) = 


vin) = 2 G(x:,y;) u(y.) , (isn) 


Xp) see. Xn I particles or gridpoints 
Yar cress Ym 


X02 Xn { (coarser) uniform grids 


Yio: Ym meshsize h 


p-erder interpelation 
F(xi) = Zwix #(X,) + 0 (pny ¢) 
. £ (y,) = >. Wye ¥( rf) F O(td)?F ) 
v(x;) & 2 Wy (xi, Ye) u (y;) 
= a G (x: Ye) U (Y, ) UlY))=4. Dye u(y] 
ae pa Wik G(X, \,) U (Y,) anterpelatich 
(B) interpelatien @ recursion 


®eorrection: viui)svlki) + 2G 


local 


(xs, ¥;) u(y;) 


n-particle foraes 
ne 


Calculated directly : FB vy operations 


Multiscale calculation : Ch operitions 


Ec dimension, accuracy (det.), uniformity ) 


Only the local part of the forces 
should be used/updated in local motions © 


Smocth part of the force 


— in larger scale (collective ) motions. 
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# operations: O (ws**') 


N gricpoints 
O(h*) accuracy 
cd cdimension 


Coelkecvenl << 4d 


Harmonic Kernels: ~ GdsN 


Via multigrid Poisson solver 
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Inteare - Differential Equation 


Lux) = [@G,9) uly) dy 


differential 


Multiqrid Solver 


distributive relaxation: 
Ast order 45 aS 


2nd orcer #§ -2d +8 


Solution cost x one fast transform 
(one fast evaluation of the 
discretized integral transform) 
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Elements of linear multigrid (MG) 
for solving elliptic PDEs 


Klaus Sttiben 


GMD/SCAI 
Schlo& Birlinghoven 
D-53757 St. Augustin 
Germany 


Sera 
nients ae 
Ea Sacer nae Meet, 


wy ESS 
ets 


ee es aks 
SHURE 


arse xe, oA Lee 


imethy 


itatior i ee 


MGi-1 


3F 


What is multigrid NOT? 


A particular solver 


What IS multigrid? 
A general strategy for constructing solvers 


This strategy exploits the fact that a problem 
can be approximated on different scales of 
resolution ("grids", "levels"). 


wat 


‘Compute paly, those "5 
oe eae which REQUIRE 


Features of multigrid 


(1) Provides "optimal" methods (i.e.: O(N)) 

(2) General applicability 

(3) Involves only "local processes" 
(Parallel computer architectures!) 


MG1-2 
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lassical methods. 
increasing increasing 
efficiency generality 


Sadcaicee i t ta ee 


Gau&%-Elimination 
Gauf£-Seidel 

SOR 

ADI 

Buneman 


Poisson equation: 


symmetric differencing (2nd order) 


unit square: h=I/n 


#unknowns: N= nZ 


Solution up to discretization error €= O(n?) 


method: ‘# operations — time, h=1/256 


Gauss-Elimin. | O(N”) approx. 1 day 


BIZ 


SOR O(N™~ ) loge approx. 30 min 


ADI. O(N logN) loge 4 min 


Buneman O(N logN) 15.5 sec 


MG (iterative) | O(N) loge 18.5 sec 
MG (FMG) O(N) 7.6 sec 


log ¢ = O(log N) 
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, Local refinement 
+ Local coordinates ~*~: 


"Finite differences 2... 
Yo. > Finite volumes. 8 2. on), 
"Finite elements 72... 
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However, .... 


... the development of concrete multigrid methods 
for complex problems may be non-trivial: 
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A multigrid tutorial, SLAM, Philadelphia (1987) 


An introduction to multigrid methods, 
Pure and Applied Mathematics series, 
John Wiley and Sons (1992) 


Multigrid Methods for Process Simulation, 
Computational Microelectronics 17, Selberherr (ed.), 
Springer Verlag, 1993 
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Ist: 
2nd: 
3rd: 


4th: 


- ; Co ‘oppe er. : Mountain MG Conference. = 


cme sBoae Re Sig oe wey se ere ese ele areas etnies inners ee Selene 


fo amass opebeeaed wet ‘oe Wa eects od ate be dead 1 ota eae nb eh eens oe te wi Alea eo aesawee: ee Nee 


Appl. Math. Comp. 13 (1983) 


: Appl. Math. Comp. 19 (1986) 


Lecture Notes in Pure and Appl. Math. 110, 
Marcel Dekkar (1988) 

SIAM, Philadelphia (1989) 

Communications in Applied Numerical Methods; 


Special Issue on MG, Vol 8, No 9&10, 
John Wiley & Sons, New York (1992) 


: NASA Conference Publication 3224, Parts 1&2, 


NASA, Hampton, Virginia (1993) 
2? 


aoe ope te owe awn ago 


"European MG Conference | 


te ee er 


Lecture Notes in Mathematics 960, Springer (1982) 
Lecture Notes in Mathematics 1228, Springer (1986) 
International Series of Numerical Mathematics, Vol. 98, 


Birkhauser (1991) 
International Series of Numerical Mathematics, Vol. 116, 
Birkhauser (1994) 
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AN 


Relaxation 


Coarse-grid correction 


Nested iteration 


First step: 
| | Multigrid | iteration Cycling") 


Relaxation + ‘Cbarse: eed correction 


Second step: 
Full multigr id (EMG)... oe 


shales gases Beans etch aa Jia cereus Pan 


Multigrid cycling + Nested eepadionl 
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AS 


Basic principle: 


Given: AX=b Splitting: A=A, +A 
One iteration step x —> X : 
Ay xX + Ag 2 deal 8) 


(-A} Ad) = asympt. conv. factor (spectral radius) 


p 


“yy pical PD a pieiadon: ae | Tate 


Bee one a cemekd ia bates totes det me tmed eines mae 28 


Jacobi relaxation: 
A, =D, Ag=L+U 
Gauss-Seidel relaxation: 
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For all grid points do 
Z + 2 + A + 
oe a Mies ie ici | 


Effect on error vie ue. ue 


h ‘ h i: 
ij =3[v a an a | Vit | 


m=) averaging process 


The error gets smooth very quickly! 


mb Peete 
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Error after 10 relaxations Error after 15 relaxations 
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Fourier ssa Ue ore error " (1D) 


screscery eee vane. a veer my ape simmer py ematy ew amrce ae 


doe frequencies or a 


‘ a , 
' be arrester 
to) oe BP Atel BTpig ete 


i eee 
(wavelength < 4h) ‘ 


High frequency reduction factor per sweep 
("Smoothing factor" LL ) 
Independent of h! 

Analysis by (local) Fourier analysis 
Poisson equation: GS: [L=0.5, RB: [1 =0.25 


Efficient smoothers exist for all "sufficiently elliptic" 
PDEs (or systems of PDEs) 


Typical reduction: 1 order of magnitude in 2-4 sweeps 


MG3-4 


ao 


Pern eer re 


ern ee eee nee eee enone eee = cam 


jerarchy of discretizations:,. 


th 


MG2-2 


FO 


fine grid: OQ. coarser grid: Q 


tion’. 


apres esa = = 


ansfeér.from,’ 9, 


7 


‘prolon gation. 


ote owed ron “ 


ore tees 


Restriction: 


Qn say ("full weighting") 
| ae 

h 6 mer: ae 

LZ ck 

Prolongation: 


(linear interpolation) 


h 1 2-74 
oh = "Ae 242 
1. 2.4 

MG2-3 


Coars e- o-grid ¢ corr rection 


l a 


h 


h 


Given: L? u" =f5 approximate solution: u 


Correction equation: 


h 


roy ca! with a®:s ¢2- Lu” (=residual) 


ay Ube why vB 


This process, applied iteratively, 
CANNOT converge! | 


Hoh , 
a d= = 0 => no correction! 


at 
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ier. components » 


“Coarse-grid-aliasing of Fou 


ais aomeie mat 


Assume: Fourier decomposition of error 


a Fourier frequencies (1D). a 
low frequencies. a - i “high frequencies — 
(wavelength | > Ah). a vies «(wavelength < ; 4h)". 


-Visible on coarse e grid! ‘Invisible on coarse grid! 


_ 6000-00 do 


a Ce 
en 


—o—e—e—-e—0—_0—0—0—0—_ 


High frequency Fourier components 
CANNOT be corrected by coarse grid! 


However: 
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~The (iterative) two- o-grid cycle | - | 


Dine ee emer e cate or ences: 


ieeleyanon + ene correction 


reduces high frequency reduces low frequency 
error components error components 


Vy: pre- ee V2! post-smoothing 


Convergence independent of h! 


Vy V2 small (typically 1 or 2) 


Rough convergence prediction: smoothing factor 


Analysis by (local) Fourier analysis 
However: Two-grid method is not yet efficient! 
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a 


mt, 


Notation: 


\ restriction 


(v) y smoothing steps 


© solve grid equation 4 prolongation 


Notation: Vi yy, V2) -cycle, W( Vir Vz) - cycle 


Compromise: F-cycle 


MG4-4 


VW := 1 relaxation sweep on finest grid 


V=V,+Vy, standard coarsening 


- 43 ae 


v=3: ow _ Se 
vw yaw ie 
2 vw AB. vw 


6 WwW 4 WwW 


1t ETEOr mared ed 
xT ae nua oe pe on. Sy Eas Rema 5 
lenta > OFLU Ki 


with Moreau 


Solution up to discretization error © 


O(N) loge = O(N) log N 
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2 Standard:MG for Poisson's equation. | 


Aer . : eS ne at 
senda atm calcined annette ae mare hint ters atti fee devtninmematnmetnnn, mein Crees tafe Sn sae 


Poisson equation 
“1. 
“5\-1 4 -1| U =f 
-1 h 


h 


" Relaxation: GS’ 


ing 


ae ae - full weight 


i “Oh 
linear interpolation 


Convergence factors (h=1/128, 7 grid levels) 


1,0 0.500 0.395 0.400 
Bae | 0.250 0.190 0.190 
1 


Zi 0.125} 0.117 0.120 


MGS5-2 


Rae. 


een ae een nr 


areal meen ee eae —) 


iry on MG comp 


onents 


SE RIO DC Tae 


Discretization MG straightforward if "sufficiently elliptic" 


Relaxation Gauss-Seidel type (point-, line-, blockwise); 
collective and distributive relaxations; 


[LU-iterations 


Coarsening Standard coarsening (h,2h,4h,...)7 
strategy Semi- or adaptive coarsening; 
Multi-coarsening; Variational coarsening 


Coarse-grid Analogous on all grids; 


operators Galerkin operators: Lie L ioe in 


Prolongation Linear interpolation; 
"Adaptive" interpolation (e.g- discrete eqns.) 


Restriction Local averaging (e.g- full weighting); 
Transpose of prolongation; 


"Adaptive" weighting 


Cycle-Typ V-, F-, W-cycle; "Adaptive" cycles 


Vy" V2 Typically 1 or 2, usually Vy 2 V2; 


dynamically chosen 
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:d MG'for Poisson's equation. 


“13 Residuals donot. - 
ace! ge much locally: 


Cree ee 
. 


0.487 0.809 
0.118 0.173 
0.034 0.059 


ed 


Pe er ee 
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ryt ae: “at 
Lay, coarse: grid | corre eehip on. 


: bi courect 
We. soe oe 


minimize nh RNs 


£(urTt Vv)” nl thy") = (ur nv b 
av aR Ey ny" va a ra) 


peers Ge Con onser ative i. ; “Fe te we htm ng. 


+ epee ren FO mes me pp ren ge cece Were ee Ee mane Pet mmeyiar om Sn oens BE ENE OLA NS Renee PONE EO OTE RS ~ sytl ft 
~~. (ee tD MPEGS DEES SPL SITE EE ELS APE EE ETE TS BEd a MON Fe tam a ed ER we spear a 


bO 


“Variational coarsening 


cae dan en es ee at har bee on med te tna 28 mete wba a ke SIL ilgce cite ale cicam cae: Smencicheletcetaine dha Seite Kavarna 


I 
Given: i Qn, Qy and i 


i 
Defies Sl) and LS 1, Ly 


Galerkin operator 


, H ‘ 
y symmetric => L symmetric 


oa pos. definite == ey pos. definite 


srs Sree rose asenenirameengenens sare nmenne yo ermerrawe sen mewn avmeme tet fo. serene 


‘Variational principle Gra pos. definite) _ 


ue ah _ error before ‘coarse-grid. correction 


: 

| 

| 

H 
i— s 

| vn : ae error after coarse- grid correction 
| 
{ 
i 
| 
| 


AL tL eneigy n norm Re range(Iy 7 


Wee Se sete aes ac ode meee ates ott etanemae tt 96 cermin ne wb cated tatters 


Remarks: 
Applicable to non-PDE problems 
Only applicable for linear problems 
The computation of LU is relatively expensive 


Origin: finite elements 
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Gf 


ann TI ee nace a Pe porn coe coe s etnemmep meni ae mee 


sep ene er rr een 


, Proo of of the variational. principle 


uh whs. := Pub wi) energy timer product 
= SY p 


+ := orthogonal complement 


R= range( Is) ; 


Two-level h h ee oe 
correction operator: K =I ae? ) IL 
Properties: 
h, 
K issymmetric w.r.t. <.,-> 
eR mm Ki vh =o 
vig QY mp KP yh =vh 


ro) Ke is orthogonal projector onto oe 


1°. «" is orthogonal projector onto 2 


— |e = [bevel 


|v» 7 Gus «Ky vel 


tH 


it 


min | vi 4 5 || 


See 


MG5-6 


(Cc 


The performance of a concrete multigrid method is 


not sensitive to modifications in the problem, unless 
its "nature" changes essentially. 


Some typical model problems 


. Poisson equation — eee Uyy =f 
Anisotropic equations | | #BUgye Uyy =f 
| | O<e<<l or E> 1 
Discontinuous coefficients + (auy)y - (buy)y =f 
Singular perturbations  -eAu +au,+ buy =f 
| O0<exl 
Indefinite problems | | — -Au-cu =f 
Nearly-singular problems es0 
Eigenvalue problems 
Higher order equations AAu =f 
MGG3-1 
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ETT ee penetra ROT ae oe op aa tenes tno netene Sete mene ti ncaa mgeretert an murene 


L = 9 -¢ 2(1+€) 


-1 
h 


oS . Gauss-Seidel relaxation. . 
For mall ed Sint do 
ee hf p+ eu Leia a +u 

ij i-1j i+] j ij-l° “ij+l 


ae [ bes 


Smoothing analysis (standard coarsening): 


w—-> 1 if e—-> 0 or E—> 0 


Behavior of error vo Pie uw during relaxation: 


h h h 


Vv = J ae ieee + V + V 
ij 2(+€) L i-lj it1j © “ij-l° ij+1 
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Influence of (pointwise) Gauss-Seidel 
relaxation on the error 
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(t he n\\ 
ae AN, Ci Ly, Shy 
KK \ aS AA AY \ 
T. yj a 
NY 
HN is oN ae | 
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LR aay, HA 
Error after 10 relaxations Error after 15 relaxations 
Weeg7 


direct 10n: 
be roy Riess ras i Nis ahs LEN tet ay 
aC 


ae 


ent0 nae. 
Hedin 
ect 


% Bee itl 


LULL 

TTT TTT 

TTT 
 y-line, TTT 
coarsening [TTT 
see EE 


es -E 21+) ~€ — -4E 2(14+4€) -4E 
b2| 4h 


-1 -1 


The anisotropy gets weaker on coarser grids! 


What to do for arbitrarily varying €=& (x,y)? 


e Adaptive coarsening 
Technically rather complicated (cf. AMG) 


"Multiple" coarsening 


Use more than one coarse grid on the same level, 
and combine the corresponding corrections properly. 
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i 
sects nie See 
Tots GORE 


Se 


Ea 


ip LILES) RR 


ee nglyicouplec 
yin Ss 
Ae a pithedivec = 
If 0<¢<l: Gauss-Seidel y-line relaxation 


If E>: Gauss-Seidel x-line relaxation . 


_ What to do for arbitrarily varying. €= -€ Guy) ? is 


Alternating line relaxation 


One step of x-line relaxation followed by one 
one step of y-line relaxation. 
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Influence of Gauss-Seidel y-line relaxation 
on the error 


Ml 
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EUS Gs 
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Error after 10 relaxations 


6S 


Error after 15 relaxations 
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Multigrid Monte-Carlo and 
Stochastic Coarsening 


OBJECTIVES MS. 
Energy minimization | | 

e ‘Near- minimum start OMG, AMG 
| Far start ms annealing 


° Homogenization 


Equilibrium statistics 


° O(n) per sample ms 
e OC) per sample Monte —- 
Carlo 


° Homogenization 


Dynamics a 
r Large At | | MG, AMG 
e Very large At ms annealing 


Le 
Stochastic dynamics 


ms. Monte-Carlo at each time step 


| 2 
Ee $2 aij (“5 ) 


Dd Fu 


ECU) = min E(u) 


O = ae = 4 2 orm . Cup-uj;)- f. 
Oui | iy ht at 
ies a;; ( Ui) =4 
h* <j,i7 


— E(U) = min E(u) ; 
Eta SJ lF (ux +uy)-fu | dxdy 


a 


Physical distribution | 


-BE(u) 
Pelubeb eh, pags 
| _BE(u) 


Z=Ze€ 


AM 


Observables 
E(u) | 
M(u)= Zu, Mul 


Averages 
<M7 = ve ree) Mu) 
Zo weer 


How to compute 2. 


Monte~Carlo 
we (Mc) 
Markov chains 

es, ul? | afi Tg ane 


Py (us P (ul =u) | 
= 3 Pres (0) Prag (ule ue) 


P (u) given (e.9., random ) 


Detailed balance (DB) 
Pe, (usu) _ Pe (w) 


a a 
Pr (uu) Pz (u') i 
—> lim P,(u) = p. (u) Fu 


k 
<M> = lim an > Mu") 


k~ oo n=1 


Metropolis 


current U <> one candidate Bog 
- Plu ua’) = min (Leone 
= Dekailed Balance 
Point ~by -pornt Metro Polr's 


extremely slow to produce a. New 


(ndepen den t SOM ple 


TH 


A. Simple Energy Basins — real ui 
Stow to equilibrate — h 
sand be sampled : Tu yh 


Coarse Morte - Carlo 


| ceil pote ah: « 
Hamiltonian: E (wTav") 


= E*(v") 
Multigrid cycle : 
‘on 1. Morte -Carlo passes 
each 2. ¥° cycles on coarser 
level : 3. Monte-Carlo passes 


Near equilibration and decorrelation 
in one cycle (¥ =2). 


Small work on coarser grids CG ae aed +) 


ca 


Independent sampling 


on coarser levels 


i. Large-scale statistics: 


accumulated on -ccarser levels 
in J = 2° cycle 


2. Vast domains: - 
only on coaVvser levels | 
by domain, replication 

(nstead of $= 24) 


3. Macroscopic dynamics : 


cheaply, direct ly , 


Bm 


Kamiltonian Simplification 
H(u) = H’(u) + O(v), a —> U 

4. Aprioré lottery 7 
Delete A tM proba bili ar ki p, (a) ~ sae 
Freeze A “mw probabi Ir ry 4 - P, (&) | 
peels 


Rllows coarse dynawics, 


Freezing should correspond Co Ane 


2. Aposteriori lotter 
es Ai 4 if a(@)>alv) 
Weedeat @ oa uw) = i oe *(a) — A(u’) 
Global convenient | 
Requives small A(u) foy~ smooth i=, 
Prohibits coarse dynamics $ | 
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B. Ising Spins Uz = +1 


E(u eae Ie Pla) ~ eT 
e Coarse moves at Hise flips 


e Standard blocks are improbable 2 
e Blocks Should have uniform sign . 


© Basin blocks OVY current UV 
disrupts statistical Fidelity 


Plu Ua) P(u.) (dekoiled balance). 
P(U— 4,) ~ P(4,) 


© Delete V,;(4) = dy UU; fees 
in probability p:; = qi bj 
Freeze V, (block e with j) tf undeleted . 


m> detailed balance. [swend sen- Wasa” 
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Process all interactions. i oar 
: ~ > 4 t 


Also. between blocks. 


Except . between any | | a 
two “coarse spins’, . | | | 
rf 
» Coarse Hamiltontan E = -2Vyuy, 


+ “ende pendent blocks". 


Repeat: Hamiltonians a Ef E* Ee? —_ 


On meshsizes 1, 2,4, 8, 


+ accumulating list of independents. 
m~cycle (Vv W, ae, 


> Multi- scale island. dynamics. 


East sampling, But no cndependence 
of coarse samplings . 


oe 


Detailed Ralance 


E(u) = -V() + €'(y) 


Probability of deleting V given ce 


=v (%)/T , tPe (OH) 
-y Z] ‘Pe (a) 


ieee freeze V 


p(G) =e 


Te Vic) A Vi): 


P(K>4) _ by (WH) Pe (BW _ pyl&) Per 
Bie) = Fa) Pee =O)” Sela) Perle) 
2 Pe (U) 
Pe (X) 
l¢ Va) = Viv): Pv (a= pylu)=p. 
Pls) _ pPe (Ku) + (1-9) Fe(Z>w) 
P(U>H) pp Pe (UB) + (1-p) Pe(u&) 


Pe(u) _ Pe (w) 
Pe (K) Pe (&) 


2- spin coarsening (Ss W) 
Q2-spin given al-ernative are pr: 
configuration H | #, H, P PB 


ee pepe To 
fete te 


3- spin 
conFrquy, 


a. f> “2” ab 
whe balance a = RE ln ae = 
reba! Salance a(-P,) 7” 


cmp 
coarse — (level 


deviation 


3-spin (p=.15) 


coarse ning 


2.- spin (sw) 


coarse ning 


asy% sites treated 


Stochastic Dynamics 


At time co position vector X" 


velocity vector y" 
potential energy E(x") 7 


Deterministic implicit time step ; 
-VE (earn) = pre 
ptt =m, ave ) >= Wen Me vee A 


<> x! minimizes Hy (xe) 


H, (x) = £E my, (SV + EO) 
Stochastic Time step : 


P (yt) me Ene Hn (x) 


iy flr 
Pe i 
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Results of the stochastic step : 
@ <6° "> S-VE) 
oo <(pttty > -<p"> =m,/B 
e Local components, with 


oszillation period K St 


are Fully thermalized: 
P (xml) w 7 Et) [kg 


© Large scale components , with 
oscillation period > St 
satisfy Newton law. 


Stochastic , Hamiltonian Stmpli Fieation 


Coarse- 


STOCHASTIC 
HAMILTONIAN 
SIMPLIFICATION 


configuration 
H, = Original Hamiltonian 

obtained from fine gril => complicated 
H.= Simplified Hamiltonian 

FAS: original form + polinomial terms 


He ia H, ~ Ho 


Detailed : P (H,— Hs) = e, Hem Hs 


Balance: 


current 


— P( freeze Hs -Ho ) = 1-P(H,— Hs ) 


‘Simpler : Approximate detailed balance 


$b 


Thermodynamic limit 


Statistical average of 
ON infinite system : N-» 00 


Calculation to accuracy Eo 
=> Large N = rt points 


xX ne Monte-Carlo sWeeps 
to produce a saniple 


x O(e~*) samples 
E calculation of 5 fog det A| 


= O(n?*4 e~*) [x oo] 


MG calculation in O (e-?) 
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Statistical Fields 
Probability (4) 4 ow Energy) PT 


Monte Carle a particle s by = particle 
(or- point - by - point ) simulation 


Slow large-scale samplings 


A. Small changes per sweep 


ce Few samples per con Figuiration 


Multigrie Monte Carlo: 


1. Moves on all scales 


ae Much sampling on coarse levels 


E.g., Gaussian model on creak lattice: 
O(A) instead of oi) operations per sample 


==> Homogenization 


4% 


A. Brandt 


Global and Discrete-State Optimization: 
Multiscale Annealing 


min E(¥) 


False attraction basins 


Particle by particle minimization 

trap ped in local attraction basins. 
~E(u)/T 

Simulated annealing P(u) ~e d/ 


T+ O 


trapped in large - scale attraction . basins, 


Multi-level annealing 

Collective moves, at all scales 
converges (fast, in probability) , | 

o A large-scale move is decided only after 


optimizing around it at all finer Scales: | 
@ Fast annealing at each level. 


6 Recombinations at each level. 


go 


OpLIMIZALIOVE — min Flv) 


a Fast convergence near optimum 


» Global optimization 
Escape false attraction basins 


E(u) 


multi-scale attraction basins 
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Malte - level 
Stochast (Co 
Optimization (minimization) 


Large basins & large -scale 
oO spatial features 


Large scale changes 
decided (stochastically ) 

only after . 
Optimizing around them 
at all finer scales 


Recursion 


: 52 | 


a 


_ : 


° 


hi 
a 


o 


«a 
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-€(4) /T 


P(u) = = & 
E(u) = ¥ (u.-4,)” ae FLU; 
<t,J> b | | 


Real U,;: f.=O = Uiw~YUj 


QO, : Ising spins U. = tt 


| a ‘ 
“a | 
+ 


Oz , O3, eee spins : vectors lu. =4 


Tslauds within islands 


Hierarchy of basins 


Many Local Minima” 


Ising sping: min EG) 


Spin glasses S; = +4. BC 
Els) = -— 2. 4y 5:55 > DS. 
Ce,j> 


Relaxation — point-by-pont minimization: 
Flipping a spin if the energy decreases 

~ BE: 5) 

Artificial temperature p(s) le . 


Point - lby- point stochastic relaxation 


Annealing Increase (> gradually 


*) : 
with Dorit Ron (1935) 


Ising Mode! 
“BIS)* “£58 j 


» Ground - State ; Siz4 


bs 
ee IESE 


Se 
mS 


GG 


4 4° 
4 4 
4 4 
iis f 
rt 4 
4 4 
“ticle” 
4 a a 4 
j 4 
L, A 
“4 rae | 
4 oa 
+ +: 
4 4 
4 4 


eS 
* External Magnetic Feild : H<Q 
*Ground- State remains’ : Sted 


Sled 


Block - by-block Relaxation 


E(s) < E(S) 


==> Local optimization around Ss 


Working Area : 


Block revision precess =» revised block 
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and highly non-quadvatic OPTIMIZATION 


1. 


2. 


Hherarchy of change scales 


Large-scale change ts decided only after 


calculating cbs: effects (optimizing around 


it) at all finer scales, starting from 


eAcceleratioy 


® 


the finest. 


: At each scale employ stochasbicity just 


large enough to eS cape local wmintma on 


that scale, then strict ment mizabion. 


; Repeat (the finer —the more) with Loc 


. Recursion 


| @ Hyper parallel 


EScaprny lecal minima sith lawny e 
ath-action basins 


GS 


Supplementary Minimization Techniaves 


4. Lowest Common Configuration 2 Lee 


¢ between the revised 
and unrevised block 

vat repeating any level 

e within the revised 
block relaxation 


2 ewer starting B 


9 increases the probability of breaking long blocks 


3. Ado 


®» Saves work 
> relaxation may expand beyond original work Areas 


tive. relaxation 


4 Pacey, chifte 
9° half size of the blocke 


5. New ty pe of blodas : broken along weak bonds 


{oo 


EXAMPLE: A complete multilevel cycle 


level 


{x1 |O} o 
finest 
oie: 


2x2 | | 


eg | 2 


@ flipping a square of spins 
stoshastic relaxation 
o minimization relaxation 


Complexity : O(N*) ; never exceeds 0(N*) 


Linear Programming 
Transportation Problem 
Lumping 2 neighboring 

| destinations Origins ; ... | 


blocks, super- blocks 


Coarse costs 
=fine marginals 


EMG- like solver 


Total worR 
z 4 simplex (1001) step 


Role. at rela xrabvon [Kaminsky | 


AVLTI" LEVEL OPTIMIZATION TECHMIQUES. 


Zquak LONS + Opt (mization 


A-FMG : | Pav neree optimization on coarse 


| Local optimization at velaxatron 


Tutevactive design: Re-solving mostly on coarse 


—O pk canal 


Z } 
| re conkee! 
Fup Wy ch 1G Wee | ext Ve wi ZACLOV 


Quadratic Gunetia wel (= lineav~ equations] 


Constratned quadiatre : FAS . Win 5 ot 


= 
Auk <4 


Non quadratic (few extrema): FAS 


Herahlu non-quadratr mult( 
oie dat i - ! (many extrema ): level 
Discre be states Us +4 3 ; DAL, hylts 


Linear programming transportation pire blew 
< = 


<amins ky 


{03 
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Full multigrid (FMG) 


and 


nonlinear multigrid 


Klaus Stiiben 


GMD/SCAI 
SchloB Birlinghoven 
D-53757 St. Augustin 

Germany 


MGG1-1 


Problems with iterative methods: 
First guess? Stopping criterion? 


Soprommation error discretization error 


(of current approximation) ~ 


Generally, it makes no sense to iterate the 
algebraic error such that 


Jur uh] «< lfu- orl 
algebraic error discretization error 


FMG is, formally, just a combination of 
MG cycling + nested iteration 


However, «... 


jo 


MGG1-2 


Sequence of m increasingly finer grids 


Oy Quy mr Qa = Qy) 


m 


For k=1,2,...,.m do 


If k=1 then ieee denotes higher 

Solve LX uk=£* order interpolation 
else (FMG interpolation). 
k-1 


k _ yk 
u-= Ty u 
k «>0: denotes the 

).- number of MG-cycles 


K 
i < MG (ak Ps i 
endif . tobe performed. 
enddo 


Q. ie Z7': EMG interpolation : . 


0. stands for, cycles 


Order of the FMG interpolation: 
atleast p (= order of discretization) 


Standard coarsening, 2nd order discretizations: 
cubic interpolation 


MGG1-3 


10 + 


is Lye a 


“discretization ertor | ae current approsimation error 


i. was. Be cineaend tees Macey tee 


Goal: en” E,. 


Assumption: E,- rE, 


oe pap tene eet Yee Oe BON Cees SORE TE em ER LTE emTINe IFS 


For. ‘each. ky the’ abs apptoanation uk is accurate up. to 
the. level of the discretization cao i. e. e ei E,, . if 


Sedeeee eS ey). igtalie Sree leeene cece ew eeetavers 6 9 eeele veya er erecehsievereierececeteiely = Ng 1810 “B08 te ee 
. 


cae meee at Peewee ena ererseeseee® 5 Cove ccervcscesecnesee® jg- j= j= “eeeeaceee® 


wot eeeweses F ,eereeeseee 


were wre ees 6 se ee mrene® 


stands for K cycles 


MGG1-4 


jo3 


ignored: FMG interpolation 
FMG is true O(N) method 


Poisson equation: ° - 
Ve,1)- -cycle, 0 = 0.1, r=0.25 —» K=1 


Wisc * = 4/3 Wry cte =7 relax on finest erid 
Solving Poisson-like problems with FMG 


costs approximately the equivalent of 
7 relaxations on the finest grid! 


k 
L U =f should "make sense" on the coarser grids 


Otherwise, start FMG ona higher level or 
use "averaged equations" on the lower levels. 


MGG1-5 


105 


Nonlinear co Uae aa approximate solution: ub 


LP vie ab makes no sense any more 


Linearization (e.g. Newton) may be complex 


The inner + outer iterations have to be matched 


_"Direct"' ‘approach. 


No global] linearization! 
(Local) linearization only in relaxation 


==> No matching of different iterations required 
==) Same cycle structure as for linear problems 
=) Sometimes easier to converge 


Multigrid performance as for linear problems 


Easy incorporation of "special techniques" 
(adaptive multigrid, ¢ -extrapolation, .....) 


MGG2-1 
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oe cate By tue ah) = at 


> Bye + why = Mahe eC Hla) 


ee ‘i : bs 
=! i . =! pit 


vi 


VH H 
e |. isatransfer operator which may differ from a 


e Note that uo and f H are defined above. They are 
NOT the solutions and right hand side, resp., of the 
H-discretized continuous problem! 


MGG2-2 


Att 


eats VE Sh seetie es te A 


LETT NA A OIE EE EOS oe erence peemeregreen ee cere Fe meee mart pene oe rere rerenymmrmerg® Sin mene 


‘Only’ smooth ¢ quantities” (ie. errors) must be interpola- 
ited back to the ‘fine grid! ‘Do not interpolate ut ‘itself! | 


bela tne ae ee nee ee aeetete Meee meant tetee seman meme mitens Seed mts mene e: it See te ee team seer Sierra rameter cememeee een mam sm oe 


Smoothing requires a nonlinear relaxation scheme 
Extension to "real MG" as in the linear case 
© This nonlinear cycle coincides with the linear one 


if the problem is linear 


o At convergence: 


h H VH Lh 


weU® and Us 1, U 


2 Vv 2h 
Standard choice: I = straight injection 


At convergence, all coarse-grid solutions are 
identical to the finest-grid solution 


MGG2-3 


AL? 


—AU + e(x,y,U) = 0 


Assumption: c Gxy,U) > 0 


Pe ee ee 


h oh h  -h 
Gr eg ed eS ee 


V 
Picard linearization (only good if c Gay, U;, , "small"): 


», omy, By 5) be OH aj) 


_.. Newton linearization: . 


hh oe 
iJ : ee ) 


: ¢(x,7, Ue. 
AS 


) as oleFs. Uy? * C069, U 


m= Same efficiency as Poisson-solver! 


MGG2-4 
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Under certain assumptions 


a 
Ry ae’. 


Applications 


T -estimation 


F le: ~ ~ / (e-1 
or example Try Cte =) Ty, T (c-1) 


Stopping criterion: Pu? -f | < | T, II 


qT -extrapolation 


H 
t t+ <= 


H ~~ 
h ol Y et 


h H 
Adaptive multigrid (MILAT) 


MGG2-5 
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and 


Howard Hughes Medical Institute 


AWS 


In molecular dynamics, calculating the 
N-body forces of interaction 
dominates computational effort. 


v= > kAl 
bonds 


+ S> kph 


angles 


+ >> kg (i tcos(nd + 6)) 
dihedrals 


oe am RS. 7 | 


atom pairs 27 


e Standard explicit integrators require one force - 
evaluation per time-step. 

e Fastest motions necessitate small steps. 
Unfeasible number of force evaluations 
needed to simulate slow motions 

over time intervals of interest. 


Mo 


Improved efficiency has been sought by: 

e Constraining fast interactions 

SHAKE .. | | . = 
reduced Variable treatments—Rice & Brunger, Moldyn, Oren Becker 
e Dividing the forces by time scales 

multiple time-step methods—Humphries, Byrne, Watanabe, Karplus 
e Super—stable implicit integration schemes 
LI—Schlick 

e Reference systems for the fast motions. 


‘LIN and related methods—Zhang & Schlick 


Ve RD 
bonds 
+ So kpAe 
angles 
+ > kg(1+cos(n¢d + 6)) 
dihedrals 
ACE pe: Bis 
eR eee 
atom pairs 7 “J Re, 27 


tiie 4 


The Equations of Motion 
The Hamiltonian system] _ 
Soe a oe 
p = —-VV(q) 


Langevin Dynamics 


{= Mp 
p= H=VV Gg) WP 


Overdamped Langevin Dynamics 


(Brownian Dynamics) 
eS Mp 
= VV GG) pr’ 
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Traditional Molecular Dynamics 
The Hamiltonian system 

¢ = Map 

_p = -VV(q) 


Discretization: Velocity Verlet Method 


| h 
Pn+1/2 = Pn— VV (an) 
Intl = Int hM~*Dyt-1/2 


h 
Pie: = Pa+1/2—5VV(Gn41) 


Verlet is: 

e Explicit 

e Second order in time 
e Symplectic 

e Time reversible 


e}/ The gold standard in MD 


us 


Constrained Dynamics 


The Hamiltonian system with constraints. 


¢ = My | 
p = -VV(q)+9'(a" 
0 = g(q) 


Discretization: Verlet with SHAKE 


h Porty- 
Prti/2 = Pa-~5VV (qn) + 59m An 


Int. = INthM*prsi/2 
0 9(9n+1) 


h : h 
Poti = Pn+1/2—5VV nti) + 5Int1 \n+1 


A typical bond constraint between atoms z and 7: 


my (oP eae 2° 
GAD) = (72, — 2) = 0: 


Torsion dynamics, multibody treatments, and other variable re- 
duction techniques use generalized coordinates which rules out 


Verlet integration 
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Multiple Time-—Step Dynamics 


‘The split Hamiltonian system. 
q = M~*p 
p = —VVeast(q) — VVstow(@) 


Discretization: reversible-RESPA 
qo = Qn 7 
| kho-+ 
LO. = Pn — “> VVsiow (Gn) 
forz7=O:k-—1, 
_ h 
Pi+1/2 = Pi — 5V Viast (a) 
Gti = git hM~*pz 4472. 
h 
Pat = Pi4+1/2 — 5V Veast (G41) 
end 
Qnt1 — 4k 


kh 
Pntil = Pk — % V Vsiow(Gn+1) 


LIN: 
‘Langevin Dynamics with Taylor series forces. 


= Mp 
[-VV (qo) - H(a-90) —---]- YP +t 


MS. 19 
| 


Discretization: 


LIN family of methods 
[solve] the linearized Sn over bo, h], 
with initial conditions qo. = 04 Done = = Dn 


dref Joe M”  Dref 
Pref = —VV (qo) — H(Gret — (6) oe axa r 
the residual with timestep h 
(q — ref) at M~*(p — Pret) 
(p — Pref) = = VV () VV(qo) + A(Gret — qo) 
—7(P — Pref) 


AX 


_VERLET - type hoger 
dominate (ea 


HACRONOLECULAR C Blo-) DYNAMICS 


*% one-step integrator 
Cone evaluation of forces per timestep ) 


* relatively high order | 
(“feels the wrveture ) 


* +ime-reversible , simplectic 
Caqood long-time behavior ) 
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“STABLE INTEGRATION ' 


- ia biomoleanlar dynamics 


see = | Bee 


Ejnitial 


AE < 0.00\-0.003 


2) = 


| Rule of -thumb: at= Y20 
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oH AX KE Saat! Ciccotti, Bexendsen, 
J. Comp. Phys., 23, 327 C1477) ) 


Ye BOND LENGTH CONSTRAINTS - small effed 


on tong: -time nora 


; ca RNGLE CONSTRAINTS - reduce mobility 
a the system 


RM.S. FLUCTUATION, pOIHEORALS (DEGREE™) 


RESIOUE NUMBER 


Figure 7. Root-mean-square fluctuations of the dihedrals angles 
: ¢,¥ over 25 “ps”: (—) NC run; (---) LC run; (---) LAC run. 


(van Gunsteren, Harplus, Nacromolewles IS) \s2e ( 1282) 
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——  Constronts (NC) 
10 ---— Lengh Constraints (LC) 
sissacetase Length ond Angle Constraints (1. AC) 


R.M.S& FLUCTUATION, Cq ATOMS (A°) 


i¢) ie) 20 foe) 40 be @) 60 
* RESIDUE NUMBER 
8 % ‘ s s 8 
4 A a. 


Figure 3. Root-mean-square fluctuations of C, atoms over 25 
“ps”: (—) NC run; (---) LC run; (---) LAC run. 


22 Phe @(N-C,-C) 


SPECTRAL OENSITY (ARBITRARY UNITS ) 


ops") 
Figure 6. Bond-angle and dihedral angle autocorrelation and 


spectral density functions: (—) NC run; (---) LC run. (a) 
Autocorrelation function; (b) spectral density. 


1a 


TABLE 3: Combined Stiff/Soft, Internal/External, and 
Short/Long-Range Nonbonded Force Decomposition. 
Comparison of Encrgy Conservation and Associated CPU 
Times Spent in the Various Force Routines for Velocity Verlet 
(0, = m= oy = 1) and r-RESPA Using the Propagator Given 
by eq 38 and eq 4¥ . 


"AC... Bry- My. ty my Log(BE) — Ro Tarctch Trent Tionioa Troabooded 


(0.25 0.25 1 bE =3.7073 0.0022 30.7 186.0 399.7 10085.9 
0.50'0.50 1 1 1 -3.0123 0.0087 15.5 93.3 201.2__5073.4 
1.00 1.00 1 1 1 —2.2388 0.0398 7.7 © 47.9 102.8 2579.6 
2.00 2.00 1 1 1 -1.0475 0.1981 4.2 23.6 $2.5 1326.7 
3.00 0.25 1 6 2 -2.9880 0.0102 31.9 189.1 407.6. 1090.0. 
3.00 0.25 2.3 2 2.9995 0.0102 32.1 92,3 198.7 1062.9. 
3.00 0.50 1 3 2 2.7944 0.0127 15.6 93.6 201.0 1078.2 


“log(A£) and R are given by eq 48 and eg 49. Ar and dr, are in 
fentoseconds. Here r, = 6.0 A and Ar = 2.0 A for all r-RESPA cases 
considered. Torch. Teed. Ttorvoo: 89d Teoatonded are the CPU times in 
seconds spent in the stretch, bend, torsion, and nonbonded force routines, 
‘respectively. CPU time spent calculating r-RESPA neighbor lists is 
included in Teonbonded: 


Hi) 


31 371 37? 373 174 MTs 376 
» 1000 crn! , 
Figure 6. Detail of spectral density /(+) of oq 52 as a function of 
wavenumber for “Exact” (1.¢. velocity Verlet Ot = 0 25 fs) (solid line), 
velocity Verlet Ot = 0.5 fs (dashed linc), and r-RESPA using the 
propagator given by eq 38 and oq 43 with a, = I, ny) = 6, my = 2, and 
61, = 0.25 fs (dotted linc) Intensities are in arbitrary units. 
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se RE SPA (Tackerman , Beene, Hartyna, Jer a7, 1940 ¢ 1242)) 


EFFICIENCY of the method depends on 
the eeumal split af the -forces 


( Hurnphreys , Friesner, Berne, jee 98 a8, 6868S (1994) ) 
crambin - 4G res. 


V = Mis . Vo + Veo si Vib 
+, 
Vent) WARC) 


~ ~ st 
oT, N,v, Rn,T, n,n, 0,0, = Av 


A394 


Watanabe, Karplus, J. Pry. Chem. 44, 5660 (1948) 


Singlet - RESP, | 
Ve. Vit Ve * Veg + Vib 


| Double 1 RESP | 
BPTI 


“TABLE 10: Average Total Energies, Fluctuations of Total 
and Kinetic Energies, and CPU Ratios for the BPTI 
Molecule 


CPU 


(Eveast) KE vast BE wwe : 
model (kcal/mol) (kcal/mol) (kcaVmol) log(AEmms* ratio? 
Standard 
J 1134.55 0.127 15.40 — 2.08 -L.00 
2 1134.48 0.545 15.76 —1.46 0.52 
SHAKE 

1061.12 0.844 15.81 — 1.26 0.28 
Single r-RESPA 
l 1134.71 0.290 15.40 -1.73 0.28 
2 1135.21 0.498 15.38 - 1.49 0.22 
x 0.21 
; Double r-RESPA 
I 1134.41 0.562 15.20 -1.43 0.26 
2 1135.36 0.662 15.48 mpl 0.22 
ay 0.19 


“See eq 18. All values calculated over 10 ps of production run. 
* Each CPU ratio was calculated with CPU time (method)/CPU time 
(standard 1). ‘ Total energies of sr-RESPA 3 and dr-RESPA 5 are not 
conserved. 


vv. 
Moat) + Vyp( te) 
ed en idl 


0.5 7 &.° -4c 
0.9 °-2.5 -5.¢ 
inten! nb, nb 


LLIN) (ttramg, Sebi, $Comp.cher, I, 1220183) 


Langevin dynamics - trlal eneras Guckverte, 


Averages for Alanine Dipeptide, LIN (At = 30 fs) 
VS. Explicit (At=0,5fs) Langevin Trajectories 29s 


Ey 19.7 3.4 20.6 3.6 
fe | ei] | ele 
$e —116.0 A791). - 119.2 33.0 
yp | 64.6 90.2 66.9 86.2 
TON T3542 O28 "34 028 
ree? 1.446 | 030 1.447 
tec 1.520 034 1.520 034 

UCNCa 4 ; 


Once 109.3 44 109.5 4.5 
8c. on 117.0 2.8 117.0]. 2.9 


Energy is given in kcal/mol 
b Angles are measured in degrees 
€Bond lengths are given in A 


Aree 


0 
90 100 110 120 130 


On-c.-c 
06 
04 
02 
0 iS 
0 100 200 300 
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(1dxa) ¥-or 


100 150 200 250 300 


50 


150 200 250 300 
time (ps) 


100 


50 


AS 


COMPARISON OF LONG TIMESTEP APPROACHES 


[METHOD | At (i) 
[VERLET / ee eA 


| SHAKE [constrains | 20 | 
REDUCED VARIABLE CONF. SAMPL- 2 ae 


| SUBSTRUCTURING ese 
| RESPA multiple timestep on 
5.0 


separating framework | 29.0_| 


substructuring 
local modes | 16.7. 
< Arar; > sa 


Efficent Dense Hessian Computations 
in 
Molecular Minimization 
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Overview of Presentation 


e Introduction and Motivation 

e Newton-Based Minimization and the Hessian 

e Need for Hessian- Vector Product in Iterative Matrix Solve 
e Goal: Describe Product in Terms of N-body Problem 

e The Hessian Matrix Structure 

e Decompose Product into Off-Diagonal and Diagonal Parts 
e Discussion of Complete Computation 


e Ongoing Work 


Adle 


Introduction and Motivation 


3 e Diverse applications of Potential Energy Minimization 


— NMR and X-Ray structure refinement 
— Molecular modeling 


— Long time step dynamics (Schlick) 
e Newton Methods 


— Quadratic convergence near minimum 
— Large computational overhead 
* O(n?) space 
* O(n') time with direct matrix solve 
* O(n”) time per iteration within iterative ma- 
trix solve 
e Goal is to reduce overhead of Newton method 
— O(n) space 
— O(n) time per iteration within iterative matrix 
solve 


— Include all electrostatic interactions (no cutoffs) 


13% 


Newton-Based Minimization 
Coal: minimize E(z), B: Rk" > ®. 
Expand the function E(z + z*) about z* as: 


Be+ 21) = Ble) + ORD, Lg eee) ag + Olllz'll). 


i=l 2421 j=l 24025 


Set gradient of quadratic approximation to zero: 
* * Pa 2 
OB(2+2*) _ OE(2*) PER) 20, Vo... 
OZ Ozk j=l Az,0z; ° 


Express this in matrix form: 


OE(z* a E(z* ) O° E(z* ) a E(2" ) 


021 0212 02122 OZ12n a 
dE (z*) aE(z*) Ez *) fe) z 
022 ++ ie ae eae 2 = 0, 
9E(z") oni ) p(s ) ole" mbes 
Ozn O2zn21 O2zn22 Oznzn 
or 


H(z*)z = —g(z’). 
Newton Minimization Algorithm 
1. Start with initial guess for the minimum, o 
2. Compute the gradient g and Hessian H at z”*. 


3. Solve the equation Hz =~-g. 
4. Set 2* to z+ 2*. If not converged, goto step 2. 


a 


Iterative Matrix Solve 


© Preconditioned Conjugate Gradients 
e Need Fast Hessian-Vector Product 
e Goal is O(n) Product 
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4 O(n) Algorithm to Solve 
N-Body Problem: 


az Or. aware 20 DO! Age: eee Nye | (ige 
as 25, 0 Jf Aza 0 ++ Aan | | a 
0 0 a8 fine ee Ang cs 0 Qn 
Q; 
i) 
=d ; 
wo(02)| | 
©, 
where , i 
Aig = == ; 


(2; — 23)’ + (yi — y;)? + (4 — &) 


b= DY gjAij- 
j#t 


Goal: Describe Hessian Matrix-Vector 
Product in Terms of N-Body Problem 


ar 
Ox, Oz, 

OE 
Or2dz, 


eer 
Oxry,02x, 

OE 
dy10z4 


ee 
dy20z1 


VE 
OynOzy 


Ee 
02,02, 


E 


Ee 
Oz (0x2 
ee 
Ox2022 


Chapt 
Ox 0z2 


@E 
Oy, 0zr2 
ee 
Oy20x2 


fake Ol 
OynOZ2 


VE 
02,022 
ae 
022022 


ee 
Ozndr2 


Hessian Matrix Structure 


Ee 


OFn0Zun 


BE 
dy zn 
VE 
Oy2U0In 


@E 
OyndEn 


Ve 
Oz,0Xn 
ee 
O220En 


‘ 


@Ee 
Ozndin 


Ee 
Or20y1 


er 
Orndyy 


&E 
dydy 
ee 
dy2dy1 


ae 
Oz, 0y2 

Ee 
Or20y2 


GE 
Orn Oy2 


@E 
Oy Oy2 
ee 
Oy20y2 


Ee 
OznOy2 


hac 
Ox 0yn 

Pr 
Oxr20yn 


Ee 
OrnOyn 


Ge 
Oy dyn 
@E 
Oy20yn 


(ka Ot 
r) YnOyn 


PE 
O210yn 
@E 
0220yn 


. 
. 


&Fe 
Ozndyn 


rH. 


@E 
dy0z, 
&E 
Oy2d21 


PE 
02,022 


@Ee 
Oy dz2 
fika OF 
Oy20z2 


ee 
Ozn G22 


er 


O2ndZn 
Ee 
Oy Orn 


@E 
Oy20zn 


@E 
7) Yndin 


aE 
O2z102n 

er 
022020 


Se ee | ls 


Off-Diagonal Components of Product 


Consider off-diagonal parts of upper left subproduct: 


ee .., @E \[zg 
0 021022 Oz, 0Zn 71 
02202} 0 Ox20En 2 
VE ee... Se 
OrnOz, OLnOx2 0 Ln 
Let: 
1 


g=-— —— rT 
V(2i — 23)’ + (ui - yj) + (4% — 2) 


We can rewrite the off-diagonal part of this matrix as: 


nz - aD O sda. os Mas nti2 
Oi? Qoeecee: 0 Ag, O ++ Aan || 2%2z, 
2 : : ; t2 
Ohi Oy se ne Nia Naps 0 Ene 
Let: 
Qs; = QiXi 
A = 0 A Li — Lj 
ij ~. ig 3 
“7 ((vi — 25)" + (Yi - yi) +(a- %)°) 
Rewrite product as: 
a 
zy ee halter 0 Aria ail Naya Qa, 
0 92755 wien! 0 Ago, 0 or Nap n ns 
0 0 eas ne Azna Arno ap 0 Qan 


yuu 


where 


4 O(n) Algorithm to Solve 
N-Body Problem: 


nz 0 aa 320 O Ayo s+ Ain] | % 
0 @7 0 Ao, O +++ Aan | | 9 
0 Qo ase Ine An An,2 oe) dn 
D; 
fs] 
~ siog(a2)| | 
Pn 
i 
Mos 
2; + (yi — Ys) + (Hi) 


©; = Nie. 
i yo i,j 


1S 


Diagonal Components of Product 


_ Diagonal part of upper left sub-product: 


OE nary = 
02102 0 0 val 
BPE ds Fe 
0 022029 0 v2 
@E = 
0 0 OrnOIn Ln 


is equivalent to: 


a oa 
oS aAts 0 iy 0 
qa 0 all ha 2 
0 0 ° < . , : 
oe 0 0 + gee D gAng 
jen 
which can be rewritten as: 
— 0 
M2157 32 0 ra 0 O Ayo «+: Ain |} % 
0 qnBrzen3 he 0 Ao : a Aan q2 
0 0 sees QnEn zz Ant An,2 58. 0 dn 
or 
Ee? 
M1215 7,2 0 0 con 
2 
0 @Boge 0 b, 
: 4 Py 
0 0 rite InEnge2 Dn 
Vey 


Chea g? 
Oy, Ory 


0 nz 
0 0 
qn 0 


an& 0 
nz 
0 
an 0 
Ax 2 Agi y 
Azan 
Aza 0 
Ax 2 Ax, n 
Azan 
None 0 
Ary, 2 Ain 
Azan 
Asan 0 
0 0 
qF2d? 
022032 0 
dnEnd? 
0 OL, 0En 


qnEnd? 
Ozn0In 


0 0 
ne 0 

0 Wns 

0 0 

7] 

55 0 

0 Qn Yn 

0 0 
q2 2 0 

0 In 


Ayn Ayn. 
0 Aya 
Ayes 0 
Nyaa Ay, 2 
10,0? 
oz 0y 
0 
0 
ny, 2? 
Oy Oyy 
0 
0 
aie? 
dzdy 
0 
0 


0 
92,8? 
Oz20y2 


9292.8? 
O220y2 


0 


az 0 
0 ne 
0 0 


a 

0 W2 ys 
0 0 
nz 0 

0 ne 


qn9, 0? 
Orndyn 


0 

0 
An ip a 
Dyn dyn 


0 
0 


an Uy a? 
Ozn0yn 
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Complete Computation 


2728" 
Oy2022 


GnUn 
gz 
q2z2 


Gnzn 


QnEnd? 


Orndin 


8 |@Q 


a | 


— 
rQ 
—— 


<2 |Q 


QD 
& |Q) 


| 


—_— 


—s OO 
rR ee) 


More Compactly... 


8) [a2\ | [ (al [Ay] [Az] 
a las [Az] [Ay] [Az] 
o2| [a2|| [el Hv) (el 


Pictorially... 


hes 


Ongoing Work 


¢ Implement rudimentary O(n log n) tree-code without multipole trans- 


lation operations 


e Derive multipole translation operations for Az, Ay, Az “potentials” for 


full-blown FMA 


e Incorporate code into MD efforts at Schlick lab 


e Parallel implementation 
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long- range electrostatics must often | 
beincluded 


eq. protein-DNA , 30000 atoms 


clivect method N? | 
FNIM to s. pret 3000N— 
16 times the cost ef using tutotts 
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i a ee ee a 


(fma-direct)/dir 


16 
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1117 water molecules (20A) 
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A 4Levels w/ FFT 
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idea - (partly ae Windemuiiy | 
ee ) arti ficial spli Hing of clectrontedics 


2 use ‘te init for smooth part 
e.g, 20 fsecs 


= SO% premivin for full electrostatics 


<F(t) - F(t-T)> for Oxygen 
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long- timestep methods : 
og positions, p momenta, M masses | 
FE short-range forces, G smooth forces 
—e Food + Gor) 
ideq | 
cy) sample G infrequently 
ci) solve “reduced problem’ 
involving F and G-samples 


worthwhile tf either | 
cy cost CG) >> cost (Ee) * and/or 


(vr) long timesteps w/ F ave economical 


example 1. LIN | 


example 2. ‘impulse method = - 
aka Nerlet: -1, RESPA © 


pes Bhs BCC: 

(xB y ett) airy of 
7 = M'p, 3 ox P = = = F(x). 
applied to exer P”; 

Pp oa pry At Gx) 

symplectic, veveveiile 
tf F+G.. 
@.g. Deke method 


O (at) impulse 
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Bm EBM 4 BE (CHC 
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: Fo a t F 
Get Gt oe 
At OAt Qh Aat 2At . 
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j-th level time step Cimpulse method) 
step(Tt, F)= 7 
jen then |. 
pi=p>+ 27At 6.00; 
Step (2)"At, Fas 
Step (ac-' At, Far); 
prep + 29At G0) 
else | | 3 
p:= p+ Q-**At FA); | 
~:= Ke 2° At.M"'p: 
pr=pr2r"at hr 


end iF | Plow of death Sp=0 
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% oe 


h,2h, Hh, Oe 


— 163 


oH 


reduced Vartable 
Molecular Dynamics 


Oren /M Becker , Tel-Aviv UY. 
Harvard J. 
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= () Thermodynamic proper lies. ; 
Was), Biological function He activity 
®) ae. new di a or protems, 
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| Empir eel eee 
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Numerical integration of 
Newton's ae of motion: 

omg = VV) 
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i hin A A 
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Limiled Length of Simulation 


1. Small integration time-step 
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Yo ® 3000 em! => Ti, F 10 fs 
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a N Néh~ ben ded Vv teractions 
( > Z of computation Lip e) 
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4 Constraints to shidiate high frequency 
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IY elimina in X-H ae. ; => Dace ISCO cm! 


2) SHAKE —iterofive’ 5, = nia 
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3. Reduce ie of DOF (oprovinatioat 
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propagate Langevine EO”): 
Mm. =-VV (Gn) — Hee + RED 
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propagate along normal ~ oye 
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(i) Rignc body Swmilat von. 
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ee a a cheep description 
of fester Melons: 


Varalle Reduction Techniques ; 
Generate OFT of mechanical systems 
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| Reduced Variable. NNplecular D NlAMICS 
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Bodies —  Euler-Newton EW" ee 
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ass Ca. L foun Ba clerio rho dopsia . 


Type of run R.m.s. fluctuation (degrecs) 
MD. 7 ATS 

One flex body 8.05 

Three rigid babe nite oo | 1.5 

Flex (8) flex (4) flex (8) — 10.79 

Flex (2) particle flex @) 12.01 
-Blex (8) particle flex (8) 18.18 
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Open Questio ns 


a partitionin . stra te es. 
— note genertion £ Update : 


bs accuracy r ve. efficiency. 
| re sont applicatious 


ie 4-helicel bundle . (rep) ok 
g- helreal bundle. (myog lob/n J 


Advanced ap plications 
(1) Enzyme hinge- benebing metion 


? . SOLVE , Res : 
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Algebraic Multigrid (AMG) 


Multi-level strategy for solving 
matrix problems AU=f 


Klaus Sttiben 


GMD/SCAI 
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D-53757 St. Augustin 
Germany 
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References on algebraic multigrid: 


Brandt, A. Algebraic multigrid theory: the symmetric case, 
Appl. Math. Comp. 19 (1986) 


Ruge, J.; Stiiben, K.: Algebraic Multigrid. In: "Multigrid Methods", 
McCormick, S. (ed.), Frontiers in Applied Mathematics, Vol. 5, SIAM, 


Philadelphia (1987) AMG-1 
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Ba need a 


if G eee es an | 
This ealntion is ‘the basis both for 
coarsening and interpolating 


Remember (variational principle!): 


Smooth error has to be "close to" 2 := range(interpol) 


Note that, generally, (algebraically) smooth error is NOT necessarily 
(geometrically) smooth and vice versa! AMG-3 
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“Multilevel method aby recursive ¢ application si 


; *) This construction is part of an AMG algorithm! 


AMG-4 
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Requested: splitting Q=FUC 
(C-points: coarser "grid") 


Low to interpolate smooth error e 
at F-points from those at C-points? 


to be “eliminated! \ 


has 


vr we ar Ope Se Tare ‘ 


AMG-5 
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Choose C such that forall ieF: N; CC 
i.e., F-points are totally decoupled! © 


Two-level method: 


holds’ exa 


Isa direct solver, i.e. residual=0 after one cycle! 
(After relaxation of F-points: ehe Biz range(I,,) ) 


However, recursive application extremely inefficient: 
Drastic fill-in + slow coarsening! 


AMG-6 
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jractical:coarse 


General requirements: 


heat lar a ieee tade er ob Celia Sra ee ae etgewentye 


dike x y 


“either. C- -points used i in interpolation” z 
2 or ‘strongly connected to: ‘such points . * ike m) | 


i) The # of. Cc: -points should be as. small. as possible * 


‘hese are the basic objectives which allow for an automatic 
and purely algebraic coarsening process! 


(1) is responsible for good interpolation of smooth error 
(2) is important for low-cost cycles 


AMG-7 
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* BOXMG (J. Dendy) was used in 
its most robust form (i.e. alterna- 
ting line relaxation, Galerkin,...) 


AMG-8 
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Multiscale Methods in 
Molecular Dynamics 


Ros 


MULTI - SCALE (MS) 
RESEARCH 
METHODOLOGY 


Problem size (+ atoms) = no | 
Objective : Total work = O(n) 


Research stages ° 
increasingly complicated problem 
But at each stage n—_~ co 


At each stage add. one complication , 
insisting on still obtaining O(n) 
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DIMENSIONS 
1D. -2D3 SD 


TOPOLOGIES 


Many non - bonded 
e Identical atoms 
e Identical small molecules 


° Several species 


Bonded chains : 
~¢ Identical atoms, stretched 
® Helics 


© Identicol amino acids 
Chains + solvent 


General 


OBJECTIVES MS 


Energy minimization 
° Near- minimum start MG, AMG 
° Far start : ms annealing 


@ Homogeni zation 


Eauilibrium statistics 


) O(n) per sample : ms 
© O(1) per sample Monte- — 
Carlo 


e Homogenization 


Dynamics | 
© Large At MG, AMG 
® Very large At ms annealing 


Lo 
Stochastic dynamics 


Ms Monte-Carlo at each time step 


Loy 


POTENTIALS 
Bonds 


e Harmonic (1) constant coupling 
; | (2) strongly variable 
° “Bond length 37 
e Length + angle 


e Length + angle + torsion 


Non-bond, local 


e Van der waals 
° Hydrogen bonding 
Non - local _ 
e Electrostatics : (1) constant dielectric 


(2) variable dielectric 


Combinations 


Stechastic Dynemics 


At time - 3 position vector X" 


velocity vector y" 
pctential energy E(x") . 


Deterministic implicit time step 
Fe) Sant! , 


nef SVic nt} a nt XX 
Py Si eh BV og Mc Bee = 


ea x minimizes Hy er) 


H, (xt!) 2 £E m, (Ve) + EO") 
Stochastic time step : 


PAX) ~ eran xr) 


ae 
Pe et 
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Results of the stochastic step: 
a: ph eee 

0 <(prty>-<pr> =m/A 
o Local components, with — 


oscillation period gz et 


are fully thermalized: 
n+! 
Poe) Re a )/kaT 


© Large scale components , with 
oscillation period #7 &t 
satisty Newton law. 


° Easy multiscale equilibration 


Li ( 


Relaxation 


{mintmez ation 


° Atom by atom GC ey 


e Ov ly ad forces es 


includeng local part of global forces 


@ For each atom — use éks 
nakuvral (internal } coovainates 


DAD 


> 42 
E(w=+ aij (wees ) 
<i,j?7 


=2 Poy. 


E(U) = min E(u) 


_ oF ie ; _ “ 
O= Ou, ey en AS er u; ) 
yee ie Ae Nee 
¥ Zip aig (Ui U;) $i 


<> £(U)=min E(u) 
E(u)= jf [$ (ur +uy)-fu] dxdy 


Algebraic Relaxation Theory 
Ax- b (2) CEQ. error e@2X-X= (2) 


Oma --* © kom 
eociduale yz Ae = b-AX . 
Normalized residuals Vv. = - ¥; /Ial = a.e /las 


Kaczmarz Relosstion: GS for AR'y= 
ne xe K-ar (r, / lal) = po). i 


Theorem. slow comergence < I7! «lel 
le \ is reduced pev- sweep at least by (AP/ (rox = laaf j 


[a.a;"] 


i oo ical converging errors ave special 
sa cam be oo by a ea - 
system (coarser grid). i +, bY Te oni 


Discrete h-elli ic ec S- clow convergence 
=> | re*| «lbllel & eh smooth 
| on Scale " 


ou4 


Slow to equilibrate a. 
and be ‘sampled Inv" 


Coarse Monte -Carlo . .- 


WJ~herh . 
Hamiltonian: | E (wsIyv') 


eet" 
Multigrid cycle: 
on - As Monte-Carlo passes 
. each 2. ¥ cycles on coarser 
| level 3. Monte-Carlo passes 


Near equilibration and deecpre lation 
in one cycle (y=2). 


Small work on coarser grids G tLetee) 


AUS 


Coarse -level variables 


@ Gridpotnts 
@ Subset of atoms. — 


ole 


Tnterpolation. 


@ of displacements | (smooth ) 
o Reflects lowest oa modes” 
amu mene el lian ae oer 
hepresenbabion of all these medes 


e Strong couplings ave Vok consbant 


@ Nou linear relabvons re tacned 


(to Che Sa vest needed) 


@ Iv tevpolatvon order 


© Dertvakcon of interpolation 
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electrostatics 
Now linear ty retained : 
eee park of the forces 


\s ae Ce geal “par ) 


(ond residne 
pay by cu pote (and hange ) 


cu tl coav $k ~ level cy namics 
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STOCHASTIC 
HAMILTONIAN 
SIMPLIFICATION 


configuration 
f, = Original Haniltonian 
obtained from fine grid. => complicated 


{ee Simplified Hamiltonian 


FAS: original form + polinomial terms 
H, > H, ) Ye ~H, 


Detailed  P(H,—-Hs) =e" 


Balance: 


current 


P(fecve Hels) 21S P (Hse ae) 


Simpler : Approximote detoiled balance 


Ap ostertowt acce 2 bance ch eck 
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(2 9D choim 
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couplim 4 Cae 
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| | : interaction, 


(3) 3D chaims 
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torsiom- rhs meighbor ; 


PROPERTIES OF-COARSE-LEVEL 
REPRESENTATION 


Coatse atoms: more freedom 


Bomd com straimts: 


; a X~ Caarse Poimts 
2D i\ eos fime-omey Poimt 


sp FR ‘ 7 - Extermal fiegol 
| | SS perpemolicular 
a to pflame of atoms 


NON~ LINEARI TIES 


AA 


Ty pichd reduction of stremgth of imteraction 


Om a course Ceved, in 2D. 
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Effect of BS om comvergemce rates 
3D gemeric chaim 
Emergy terms: 


K,( T-5) fi bomdimy 


* K,(e-@,) - amy Le. im teractiom 
Kcp-9) - torsiom 

N=26.. 

on 

Rate 
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-INTERPOLATIONS 
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LOCAL SETS 


For devivimg imter polation ; 
: Coefficiemt 
a: Mimimaze emevgy Localéy: 


Qocak yelaxatioms 


direct solution. 


Locat ea is Qa trum cated Chaim 


- Igmore okt-of- set limks 


Facil: aes accurate Smaath mdtigns 


Sige of locate: Set 
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